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Abstract 

Electron spin decoherence caused by elastic spin-phonon processes is investigated comprehen- 
sively in a zero-dimensional environment. Specifically, a theoretical treatment is developed for 
the processes associated with the fluctuations in the phonon potential as well as in the electron 
procession frequency through the spin-orbit and hyperfine interactions in the semiconductor quan- 
tum dots. The analysis identifies the conditions (magnetic field, temperature, etc.) in which 
the elastic spin-phonon processes can dominate over the inelastic counterparts with the electron 
spin-flip transitions. Particularly, the calculation results illustrate the potential significance of an 
elastic decoherence mechanism originating from the intervalley transitions in semiconductor quan- 
tum dots with multiple equivalent energy minima (e.g., the X valleys in SiGe). The role of lattice 
anharmonicity and phonon decay in spin relaxation is also examined along with that of the local 
effective field fluctuations caused by the stochastic electronic transitions between the orbital states. 
Numerical estimations are provided for typical GaAs and Si-based quantum dots. 

PACS numbers: 72.20.Ht,85.60.Dw,42.65.Pc,78.66.-w 



I. INTRODUCTION 



Electron spin relaxation in semiconductor quantum dots (QDs) has received much at- 
tention as they are a natural candidate for the qubit in quantum computing (QC) (see 
Refs. Q and 3 for a brief review). Particularly, a significant effort has been devoted to 
investigate the spin transition probabilities associated with the spin-flip interactions, for 
example, the so-called longitudinal spin relaxation (i.e., energy relaxation) induced via 
the inelastic process of phonon absorption/emission by an electron spin localized in a QD. 
Although generally an importance source of decoherence, the spin-flip events can sometimes 
be efficiently suppressed by an appropriate choice of external conditions including the mag- 
netic field strength or direction. In this case, alternative mechanisms (such as the elastic 
spin-phonon scattering that does not assume energy exchange between the electron Zeeman 
energy and thermal bath) can be primarily responsible for decoherence. So far, however, 
the elastic processes have received relatively little attention. This disparity may be due, 
in part, to the mathematical simplicity in the treatment of the former case (e.g., calcu- 
lations based on the Fermi's golden rule), while the elastic processes of spin decoherence 
need more sophistical approaches. At the same time, recent proliferation of the literature 
on electron spin relaxation has exacerbated the confusion on this complicated subject mat- 
ter. Frequently, an individual study addresses only a specific aspect of spin decoherence or 
dephasing without elucidating its place in the general theory or the interrelationships with 
other relevant works. Consequently, one faces difficulties each time in determining which 
existing mechanisms, processes, or theories are relevant for spin relaxation under the given 
external conditions and specific QD properties. 

The purpose of the present study is to systematically identify and analyze the crucial 
transversal (T 2 ) spin relaxation processes in the context of qubit dephasing in the electron 
spin based QC. Specifically, the role of spin-phonon interaction on the elastic decoherence 
processes is examined theoretically from the perspective of the fluctuations in the effective 
magnetic field. The calculation identifies the conditions in which the elastic spin-lattice 
processes can dominate over the inelastic counterparts. The analysis is applicable to the 
electron spin in both group IV and III-V QDs. Note that the related issues from the 
spintronics point of view were surveyed in Refs. [] and 4. 

The rest of this paper is organized as follows. Section II provides a systematic classifica- 
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tion of most, if not all, of the processes identified so far that lead to a type of electron spin 
relaxation (i.e., dephasing, decoherence or energy relaxation). Then, the discussion focuses 
on those caused by the spin-phonon interaction and describes different manifestations of 
the spin-phonon interaction in a manner consistent to both elastic and inelastic processes. 
Subsequently, the elastic decoherence processes that warrant further analysis are recognized 
and their potential significance outlined in the context of QC. Section III summarizes a 
theoretical framework necessary for the mathematical description of electron spin relaxation 
in terms of quantum kinetic equation. 

Each of the following three sections describes a particular elastic process that can be 
dominant under a certain condition. In Sec. IV, we analyze the effects of spin precession 
fluctuations due to the irregular phonon phase disturbances originating from phonon relax- 
ation via the lattice anharmonicity or imperfections. Such a situation applies not only to 
real (thermally activated) phonons but also to zero-point phonons that the classical picture 
cannot address. Section V considers spin phase diffusion due to the modulations in the 
longitudinal g-factor and hyperfine interaction through the phonon-assisted transitions be- 
tween the lowest electronic states (without spin-flip). Section VI concerns electron spin in 
a semiconductor QD with multiple equivalent minima (e.g., Sii^Ge^). As the degeneracy 
is lifted due to structural inversion asymmetry, phonon-mediated transitions between the 
"valley-split" states can be allowed when the local crystal imperfections are present. This 
process can be very significant since the energy separation of the valley-split states is gen- 
erally comparable to the thermal energy in the QC operation. Then, the discussion on the 
comparative importance of different processes in the parameter space is provided in Sec. VII, 
followed by a brief conclusion at the end. 

II. SPIN RELAXATION PROCESSES AND MECHANISMS 

For a systematic analysis, Table I shows simple classification of the relevant spin relax- 
ation processes identified so far. The term "dephasing" (with the characteristic time T 2 *) 
concerns phase relaxation of a QC system consisting of multiple qubits, while " decoherence" 
(T2) is for a single qubit (i.e., single electron spin in a QD). Hence, decoherence contributes 
to dephasing (T2 > T 2 *). The "random local fields" category denotes the processes that lead 
to phase diffusion of the multi-qubit system without causing decoherence in the individual 
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qubits. The decoherence or transversal relaxation has contributions from both the elastic 
and inelastic processes. The latter involves spin flip transitions leading to the energy (or 
longitudinal) relaxation of the spin states (i.e., Ti). Of the spin decoherence processes, those 
marked by the # sign in the table are induced by the spin-phonon interaction. The listed T 2 
processes correspond to the uncontrollable interactions with the environment that assumes 
to be in a thermodynamic equilibrium at temperature T. Inaccuracies in the parameters of 
the control pulses during QC operations as well as uncontrollable perturbations introduced 
with these pulses are beyond the current description. Table II summarizes various manifes- 
tations of the spin-phonon mechanism in a manner consistent to both elastic and inelastic 
processes. A brief description on each category is given below. 



A. Dephasing in random fields 

The fastest relaxation in a quantum computer consisting of a large number of qubits 
[N ~ 1000 (see Ref. [2])] is the dephasing process associated with transversal magnetization 
loss with the characteristic time T 2 *. Clearly, the electrons experience a dispersion of spin 
precession in the presence of a random effective magnetic field. Even when the magnetic 
field can be considered a constant spatially and temporally for a given qubit within the 
time scale of QC operation, its random variation across the ensemble of qubits can induce 
dephasing of the system. Although some of its effect may be potentially mitigated by QC 
algorithms,- 1 ^ the absence of this additional complexity provides an obvious advantage. 

One reason for the dispersion in the precession frequency is the (/-tensor variation over 
the N qubits. Consequently, even in a very pure crystal, the unavoidable fluctuation in 
the QD sizes results in the g-factor dispersion Ag, which subsequently induces dephasing 
with the rate of T* t _1 = Ag^sB (where /is is the Bohr magneton and B the external 
magnetic field). This mechanism of spin phase relaxation has not yet been investigated in 
the zero-dimensional (0D) structures. A related case of 2D hole localization in a fluctuating 
potential was studied in Refs. and Q. In addition, the randomness of spin-orbital coupling 
in a quantum well due to the fluctuation of dopant concentration was considered in Ref. [8. 

Similarly, electron interaction with nuclear spins can induce the spin phase diffusion due 
to the random distribution of the nuclear spin states and their spatial locations.- 1 ^^- In 
principle, replacement of the magnetic isotopes by nonmagnetic counterparts (/ = 0) can 
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remove or suppress this source of relaxation.— However, it may not be a practical option 
as a significant level of purification is required to be effective; the dephasing rate depends 
only weakly on the nuclear spin concentration rij (i.e., oc y^nj). As indicated in Table I, 
the random fields of magnetic impurities also contribute to the dephasing process. The 
last one listed in this category, the magneto-dipole interaction between the qubits, plays 
a specific role since it can be incorporated (thus, eliminated) to the QC algorithm under 
certain conditions.-^^ 

B. Inelastic processes 

Among the inelastic processes in Table I, the direct spin-flip via one phonon emis- 
sion/absorption (see Fig. 1) is the most studied due to its simplicity in concept as well 
as the importance at low temperatures.— J£iILI£iI2i22i The rate of this process is proportional 
to the phonon density of states g(u> q )[~ uj%\ at the resonant phonon energy hu q = gfisB 
(= hjj) as well as the square of the spin-phonon matrix element (~ u q B 2 ). As a result, a 
strong magnetic field dependence (~ B 5 ) is expected for the direct spin-flip. Thus, at low 
magnetic fields, the effect of this process can be significantly reduced. When /c#T < Hlj, 
the influence of the temperature can be ignored as the phonon emission determines the 
relaxation rate. 

The Orbach process^ 1 ^ considers the spin-flip events through successive transitions be- 
tween the ground and excited orbital states separated by an energy <5 . For example, an 
electron first makes a spin-flip transition to the excited state with opposite spin by ab- 
sorbing a phonon with energy <5 T ^> followed by the return to the ground state through 
emission of a phonon with energy 5q- An alternative channel is to make a transition to the 
excited state without spin flip (i.e., So) and, then, return to the ground state with opposite 
spin (i.e., So T Jw)- Thus, the net effect is the electron spin flip with the energy of ftw (i.e., 
the Zeeman splitting) transferred to the thermal bath. Under certain conditions, the corre- 
sponding matrix elements can exceed those for the direct spin-flip transitions between the 
Kramers doublet.— 3 It should be stressed that the Orbach process assumes two one-phonon 
transitions in tandem. Hence, the first transition with the energy threshold of So =F hu or So 
provides the primary temperature dependence along with a weak influence of the magnetic 
field. 
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The two-phonon process concerns the longitudinal spin relaxation via two-phonon ab- 
sorption, emission or inelastic phonon scattering (Raman process). It can be further divided 
into three groups; the direct spin— two-phonon transitions due to the nonlinearity of the 
spin-phonon Hamiltonian^ [Fig. 2(a)], the virtual transitions mediated by one-phonon in- 
teractions through the excited states^^ 7 - [Fig 2(b)], and the anharmonic Raman process^ 1 ^ 
[Fig. 2(c)]. In the case of virtual transitions, the integration over phonons includes the 
poles at the phonon energy t\uj q ~ 5o if 5q < ksQ, where 6 is the Debye temperature. 
The treatment of this singularity (the so-called phonon resonant fluorescence) beyond the 
conventional perturbation theory^ involves essentially the finite lifetime tr of the excited 
states that exactly reproduces the relaxation rate of the Orbach process.— 1 ^ This coinci- 
dence/similarity is lifted (i.e. poles can be ignored as compared with Orbach process) when 
the finite phonon lifetime r g is considered since r^ 1 <C r q 1 . As the processes involving the 
inelastic phonon scattering show a strong dependence on temperature (T 7 — T 11 ),— ^ it has 
been mostly ignored in the consideration of QC that operates at very low temperatures. 
However, a recent analysis^ suggests its potential significance in the low magnetic field 
regime due to the weak dependence on the field strength (-B 2 ). 

The process of directly transferring the electron Zeeman energy to the nuclear dipole- 
dipole reservoir is another possible channel for longitudinal relaxation.— This, however, can 
be effective only in the external magnetic fields comparable or weaker than the local nu- 
clear fields, which are generally sufficiently small. In addition to those listed in Table I, 
the simultaneous flip-flop of several spins generally entails energy exchange with the phonon 
or spin-spin reservoir. Much less attention has been paid to such processes in the litera- 
ture .-^^-^ 

C. Elastic processes 

The elastic processes can be imagined as a result of uncontrollable fluctuations of Zeeman 
frequency or its equivalent, i.e., the fluctuations of longitudinal magnetic field that conserve 
the electron spin projection on the quantization axis. As summarized in Table I, there are at 
least four different manifestations of this process. Of these, the first one listed ("anharmonic 
vibration") represents the case that the irregular phonon phase disturbances originating 
from the phonon relaxation induces a net effect of electron spin phase decay via the spin- 
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phonon interaction. If the phonons are treated as harmonic oscillators without damping, 
the longitudinal component of the "effective" interaction field takes the form of a harmonic 
perturbation and does not cause alteration in the electron spin phase.— [One exception may 
be the small spin phase decoherence (~ 10 -9 ) that can be acquired under the additional 
assumption of identical phase for all the harmonic phonon modes at the initial instant.—] 
However, a different situation can be realized when the phonon oscillation is interrupted by 
a series of random phase disturbances as schematically illustrated in Fig. 3. The reason of 
such phonon phase steps can be lattice anharmonicity, phonon scattering at the impurities 
or lattice defects, etc. These irregular phonon phase perturbations affect the electron spin 
precession resulting in its phase relaxation.— 

The fluctuations in the local effective magnetic field (and, thus, spin decoherence) can 
also result from the transitions between the electronic (orbital) states with different Zeeman 
frequencies. As the electron experiences the stochastic transitions between these orbital 
states without spin flip, electron spin phase can change due to the finite and random lifetime 
in the states with different precession frequencies (Fig. 4).- 5 It should be emphasized that 
although the transitions are mediated by phonons, the energy exchange between the Zeeman 
and phonon reservoirs does not take place in contrast to the Orbach process; i.e., this case 
undoubtedly belongs to the elastic processes. 

Electron spin can accumulate the phase shift through random thermal changes in the nu- 
clear spin distribution as well. Although a quantum state of individual nuclear spin generally 
has a relatively long lifetime, a large ensemble of nuclear spins (or spin flip transitions) pro- 
duces a rapidly evolving effective field distribution. Its significance can be gauged in terms 
of electron resonant frequency diffusion (i.e., the so-called spectral diffusion) in the spin-echo 
experiments.— i2£i2L2&2£ A peculiar property of this process is the cubic exponential decay 
of spin-echo signal ~ exp [~t 3 /T'ij], where the memory time Tm can be calculated in the 
model of uncorrelated nuclear spin flips 38 or in the model taking into account the flip-flop 
processes for nuclear spin pairs. 39 Note that the spectral diffusion process is probably irrel- 
evant to QC where only the initial phase decay of ~ 10~ 4 — 10 -6 matters for fault tolerant 
operation. In the corresponding short time frame, the processes with a linear exponential 
decay prevail as demonstrated in Ref. [12|. Decoherence can also be attributed to isolated 
(i.e. non-interacting) nuclear spins.— This is because each nuclear spin precesses around the 
electron hyperfine field, which is inhomogeneous over the QD volume. Such uncorrelated 



7 



nuclear spin precession, in turn, leads to the change in the strength of nuclear hyperfine 
interaction field affecting the electron spin and its phase. The efficiency of this process, 
however, is much weaker than that of spectral diffusion when the inter-nuclear interaction 
is taken into account. 

The elastic processes in Table I include the one mediated by fluctuations of Berry ge- 
ometric phase.— The efficiency of this process compared to the others is not clearly 
understood yet. In addition, an elastic phonon scattering on a group of spins can cause 
them out of phase.— ^ However, this process is less relevant for single qubit decoherence. 



D. Mechanisms of spin-phonon interaction 

Table II summarizes the microscopic mechanisms of spin-phonon interaction that give 
rise to the elastic and inelastic spin relaxation processes (marked by # in Table I). Despite 
the great multiplicity of the mechanisms considered in the literature, they can be classified 
into three groups following the global characteristics; i.e., the mechanisms resulting from 
the (i) spin-orbit, (ii) hyperfine, and (iii) spin-spin interactions. In turn, these interactions 
manifest themselves in a number of different ways, all of which can contribute to spin phase 
and longitudinal relaxation. 

Taking into account that interactions (i), (ii) and (iii) mentioned above possess, generally 
speaking, non-zero matrix elements between the states appurtenant to the Zeeman doublet, 
their modulation by phonon vibrations can be expressed as a Hamiltonian realizing the spin- 
phonon interaction (listed as "Direct" in Table II). Mathematically, if a phonon mediates 
an electronic potential V p h(r), the spin-orbit interaction is changed by the term H so - p h = 
2m 2 c 2 [VV^fe x p]s ("SOI modulation" mechanism). At the same time, the shift ARj of an 
atomic position from the site due to the lattice vibration (i.e., phonons) also changes 
the Hamiltonian for the hyperfine interaction Hhf(Rj) = X]j s ^/( r — by Hhf~ p h = 
H h f(Hj+ARj)—H h f(Hj) ("HFI modulation" mechanism); here, & h f(r— Rj) is the hyperfine 
interaction tensor. The case of spin-orbit interaction adapted to the QDs was considered 



in Ref . 
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Direct phonon fluctuations of hyperfine interaction were studied for a donor in 
silicon,— F-centers in alkali halides, 24 and adapted to QDs.— 

In a similar manner, one can introduce the spin-phonon interaction mediated by the 
spin-spin interactions H ss _ p h. In this case, one can further separate the dipole-dipole and 
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exchange interactions. Obviously, these mechanisms depend on the inter-spin distances. 
Phonon modulation of dipole-dipole interaction ("Waller" mechanism) was considered in 



Ref. 



46l ; the role of exchange interaction including the Dzialoshinski-Moria terms 4 ^ 1 ^ (that 
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can be important in crystals with no inversion symmetry 49 ) was analyzed in Ref. 

The direct mechanisms considered above are not the only manifestation of the spin- 
phonon interactions. In the strict sense, the Kramers doublets actually represent non- 
multiplicative functions of spin and space variables due to the the spin-orbit or hyperfme 
interactions (i.e., mixing of orbital and spin states; "Admixture" in Table II). This makes a 
spin-independent phonon potential capable of evoking the transitions between the Kramers 
doublet. Such a mechanism becomes apparent in the case of hole spin relaxation in a QD, 
where the spin-orbit interaction imposes a linear combination of different spin and orbital 
states for the basis of valence band.— ^ Mathematically, the mechanism constitutes only the 
non-zero matrix elements that can be symbolically transformed to the form 

'(<? Tl V rh (r) \e r> (e T| H miI \g |> (g f| H mix \e |> (e i| V Th (r) \g I) 



Eg — E e Eg — E e 



1) 



if a perturbation theory is applicable for the system with the electron energies E g and E e in 
the ground (\g)) and excited (|e)) states; f and j denote different spin states. Equation ([T]) 
assumes that the diagonal part of the Hamiltonian H mix enters into E g and E e , while the 
off-diagonal element mixes the ground and excited states with opposite spins. 

In the second-order perturbation theory, Eq. (JT]) represents the spin-phonon interaction 
as the process of virtual transition to the spin-flip excited state due to H m i X followed by a 
transition to the ground state without spin flip via the spin-independent electron-phonon 
interaction, or vice versa. In the case of electron in a QD, this mechanism can possess 
an advantage over the direct spin-phonon interaction due to the relatively strong spin cou- 
pling between the non-Kramers (i.e., orbital) states. A distinctive feature of such indirect 
spin-phonon interactions is their proportionality to the matrix elements of electron-phonon 
interaction between the ground and excited states. Consequently, these mechanisms can be 
further classified based on the types of electron-phonon interactions as well as the potentials 
that cause effective mixture of the ground and excited states. 

The mechanism caused by simultaneous manifestation of spin-orbit and spin-independent 
electron-phonon interactions was initially developed by Kronig^ and Van Vleck— for mag- 
netic ions. Adaptation of this "Kronig-Van Vleck" mechanism to the conduction electrons 



in a QD can be conveniently achieved by extracting the spin-orbit coupling in the form 
of the Rashba and Dresselhause terms (which efficiently mix the QD ground and excited 
states).—* 5 ^ Then, the phonon potential of various origins such as the deformation potential, 
piezoelectric potential,— or the ripple effect at the QD interface* 5 - 6 - gives rise to the 
effective spin-phonon coupling in the second order. The spin-orbit coupling in the valence 
band is generally much stronger; hence, holes are less attractive as a qubit.— ^ 

When H m i X is due to the hyperfine interaction, the spin flip can be considered in a 
manner similar to that with the spin-orbit interaction as discussed in Refs. 
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isomorphous case of a diluted magnetic semiconductor QD was also investigated.— It seems 
that this admixture spin-phonon interaction is more effective than the direct counterpart 
(i.e., phonon modulation of the hyperfine interaction). As to the spin-spin interactions 
(both dipole-dipole and exchange), they can also induce the mixing of spin and orbital 
states similarly; however, such a treatment has not yet been developed to the best of our 
knowledge. 

In addition to the direct and admixture mechanisms discussed above, one can also consider 
"indirect" manifestation of spin dependent interactions. This group of mechanisms accounts 
for the random fluctuations in the effective field through the phonon-mediated stochastic 
transitions between different orbital states. As the characteristics of g-tensor, hyperfine field, 
or spin-spin interactions vary with the orbital states, electrons experience uncontrollable 
changes in these parameters while undergoing spin-independent phonon scattering between 
them. Such fluctuations cause the transversal relaxation of electron spin 35 along with a 
potential contribution to the longitudinal relaxation. However, the latter (i.e., spin-flip) is 
expected to be small due to the negligibly small influence on the spectral density at typical 
Zeeman frequencies. 

Finally, there is a higher-order phonon-mediated mechanism that can affect the spin relax- 
ation. Although the phonon potential V p h(r) cannot influence spin precession directly, its in- 
corporation with the electron potential Vo(r), which contributes to the g-tensor in a different 
manner,— ^ results in a specific spin-phonon interaction H g _ p h = J2(5gij/5Vo)V p h(r)fj,BBiSj, 
where Sj are the components of electron spin operator. This mechanism was analyzed in a 
QD with respect to the Orbach spin-flip processes^ as well as in a bulk semiconductor. 61 

In summarizing the possible mechanisms of spin-phonon interaction, it must be pointed 
out that all of them make correlated contributions to the processes marked # in Table I. 
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For example, the direct spin-flip process results from the phonon-induced fluctuation of the 
effective magnetic field (with the resonant frequency fruS) that is realized by additive contri- 
butions from the spin-phonon mechanisms listed in Table II (i.e., the spin-orbit, hyperfine, 
and spin-spin interactions). Consequently, the corresponding matrix elements of individual 
mechanisms must be summed before the probability of any process is calculated. Although 



the imp ortance of interference for certain mechanisms was already discussed in Refs. 
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and 



18l it requires further theoretical scrutiny. However, the anticipated difficulty may be 



circumvented for reasonably weak magnetic fields if the spin-phonon Hamiltonian H q _ p t is 
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considered in terms of a phenomenological tensor A following the analysis of Ref. 
parameter A is designed to represent the net/total contribution of the spin-orbit mechanisms 
and can be found experimentally from the dependence of electron paramagnetic resonance 
signal on the applied stress. 

As for the contributions of different processes (Table I) to spin relaxation, they do not 
interfere and, thus, can be considered independently. As mentioned earlier, the present pa- 
per investigate dominant elastic spin relaxation processes via the spin-phonon mechanisms, 
particularly those originated from the spin-orbit and hyperfine interactions. 



III. BASIC EQUATION FOR ELECTRON SPIN EVOLUTION 

With the mechanisms responsible for the spin-phonon interaction H$l ascertained in 
Table II, the evolution of electron spin can be described according to the rubrics under "De- 
coherence (T 2 )" in Table I. In contrast to the calculation of HgL that requires individualized 
treatments for each mechanisms, most processes, elastic and inelastic, marked by # in Table 
I (with an exception of "geometric phase" process) can be approached by a general formal- 
ism based on kinetic equations. For simplicity, the derivation of these equations is restricted 
to the first Born approximation in the present consideration. As such, the two-phonon 
processes, which appear on the second order of H SL [Fig. 2(b)], are not included. 

There are several comparable approaches to describe the evolution of spin S = 1/2 
interacting with a thermal bath. For the purposes of decoherence evaluation in the QDs, it 
is desirable to apply a system of equations directly for the mean spin components (rather 
than for the density matrix components), which analyze their temporal progression on the 
Bloch sphere from an initial state to that corresponding to the thermal equilibrium of the 
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system. At the same time, it is desirable to express the relaxation coefficients in terms of the 
spectral density of the thermal bath operators that allows application of the formalism well- 
developed in the case of thermodynamic equilibrium. Finally, the equations must correctly 
account for the anisotropy of spin interaction with the thermal bath. Such equations were 
studied previously in the cases of specific symmetry properties of the correlation functions.— 
This section generalizes and adopts these equations to the problem of spin decoherence. 

It is convenient to represent the Hamiltonian of one electron system in terms of the 
effective spin operator s for the Kramers doublet {g} = {\g) , \g)} of the ground electronic 
state. The linear magnetic splitting of Kramers doublet {g} is described by an effective 
Hamiltonian in the common form H$ = g^ u pBB^s u , where the g-tensor g^ u is introduced. 
Similarly, the electron spin interaction H$l with the thermal bath (of Hamiltonian Hi) is 
expressible as a linear form of the spin operators, Hsl = S"2^s if only the states {g} are 
involved in the relaxation. Actually, Hsl also describes the transitions between different 
Kramers doublets. If these transitions from the ground {g} to excited {e} states are virtual 
(as is the case at sufficiently low temperatures kT << Sq, where Sq is the energy separation 
between the {g} and {e} states), reduction of the Hamiltonian Hsl to the basis {g} of spin 
operators s poses no specific problem. For example, the off-diagonal elements of Hsl between 
{g} and {e} states can be eliminated through the procedure of canonical transformation. 
Hence, the evolution of the spin system at low temperatures is fully described by an effective 
2x2 Hamiltonian, 

H = H S + H L + H SL . (2) 

This Hamiltonian can also be applicable to the cases other than at low temperatures. Note 
that the problems involving a group of Kramers doublets (e.g., {n} = {g},{e}) sometimes 
allow description by an effective Hamiltonian H s with {{n}\ H s \{n'}) = $n,n'{E n + 
gfflnsB^Sv) and ({n}\ H S l \{n'}) = (0) n , n /s; here, the matrix of orbital energy 8 n ^E n 
can alternatively be considered a part of Hl as well. Hence, this case leaves room for the 
Hamiltonian in the form of Eq. (|2j) even at elevated temperatures. Some examples will be 
discussed later. 

For a detailed description, we start with the kinetic equation for the spin density matrix 

a = Tr L p . (3) 
Here p obeys the Liouville equation ip = [H,p], the trace is taken over all variables of the 

12 



system except the given electron spin, and h and ks are set to 1 in this section. Assuming 
that the thermal bath is in a thermodynamic equilibrium, one can introduce its statistical 
operator 

exp (-/3H L ) 1_ 
~ Tr L exp(-(3H L y T' 



Following the approach developed in Ref . l65l with a projection operator method (see Ref.l66l). 
the kinetic equation of operator a in the Born approximation is obtained as 

^ = -i[H s , a(t)} - j* Tr L [V (t, t) , [V (t, t') , fa (t)]]dt', (5) 

where [A, B] = AB — BA. Equation (jSJ) implies renormalization of the spin-bath interaction 

v = (n<») - (n^)) s, (6) 

so that the averaging over the bath (...) = Tr £,/... results in (V) = while Hg = H s 9 ^ + 
/f2^\s. The last term in Eq. ((Sj) resembles a collision integral for the common kinetic 
equation, where the possible influence of alternating magnetic fields is ignored. In the case 
of time-independent H$, the operators on the right-hand side of Eq. (jSJ) take the forms 

V(t, = exp {-iH s (t - t')} V(t') exp {iH s (t -if)}, (7) 
V(t) = exp(iH L t)Vexp(-iH L t). (8) 

Although the integrand in Eq. (jSJ) can be calculated for arbitrary coordinates, its subse- 
quent application is not very convenient. Without the loss of generality, the z axis can be 
chosen along the quantization direction that reduces the spin Hamiltonian to 

H s = lus z . (9) 

Note that the z axis, generally speaking, does not coincide with the direction of the 
external magnetic field if the (yf-tensor reflects the low symmetry of a QD and uj = 
\J^i (di^sBi + ^uj^yj = g e ffHBB. Here are the principal values of g-tensor while 
Bi and fij are the vector components along the main axes of g-tensor.' 67 By adopting the 
specified coordinate system, Eqs. (JTj) and (jHJ) can be easily calculated with = O^ — (flffl /, 



V(t,t') = Q x (t')(s x cosm + s y sinu) + Q y (t')(s y cosw — s x smu) + Q z (t')s z , (10) 

v(t) = n x (t)s x + n y (t)s y + n z (t)s z , (ii) 
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where u = uo{t — f) and f^(t) is defined similar to Eq.Q. One can see that the collision 
integral can be expressed in terms of double-time correlation functions 

(n,(t), Q v (lTj) = Tr L {e^ H ^{t), fi„(f )} (12) 

with an evident property (0 M (t), £l v (t')) = (fi M (t - t'), n„(0)) = (^(0), Vt v {t' -t)). We con- 
sider that these correlation functions reduce to zero as soon as \t — tf \ exceeds the correlation 
time t c . Thus, two regimes of spin evolution must be distinguished (Ref. I68T ). In the first 
regime for short times t <C r c , the relaxation coefficients in Eq. (jSJ) are proportional to t. 
This leads to a quadratic-on-t deviation from the initial states of the systemthat reproduces 



the very general results obtained previously in various analyses (see Refs. |69| andl70l). Such 
a non-linearity results in only a marginal effect of spin transition from the ground state at 
sufficiently short times (Zeno effect^). Consequently, we focus on the second, opposite case 
of reasonably long time t ^> r c , when the upper limit of the integral in Eq. (jSJ) can tend to 
infinity, t —>■ oo. In this case, the relaxation coefficients are independent of time and Eq. (jHJ) 
describes the Markovian process with exponential evolution of d(t). 

To proceed further, it is convenient to introduce the Fourier transformation of the corre- 
lation functions [Eq. (JT2]) ]. 

1 f°° 

7/1 » = (fl M (r), ^(0))^ = — / (n M (r), fi„(0)> e^dr. (13) 
Considering the functional form of Eq. (jSJ), another useful relation is 

f°° (^(r), n,(0)) e^dr = 7r^(co) + iP ^ ^^du', (14) 

JO J-oo U — U 

where we use the symbolic identity L e iuJT dr = tc5(uj) +iP^; P denotes the principal value 
of the integral. Finally, let us express the density matrix a(t) in terms of mean values 
(s x ) = X, (s y ) = Y, and (s z ) = Z, 

d{t) = - [T + 4(Xs x + Ys y + Zs z )] , (15) 

so that Tr s a(t)s^ = (s M ), where a(t) is taken from Eq. ( fl5|) ; 1 is the 2x2 identity matrix 
and the trace Tr s is taken over spin states. Multiplying Eq. (jSJ) sequentially by s x , s y and 
s z from the right and taking Tr s , one can obtain kinetic equations for S= {X, Y, Z}, which 
define the spin vector evolution on the Bloch sphere in the most general form 

Y S = ^ xS-r(S-So) . (16) 
14 



In this equation, the 3x3 matrix T represents a rather complicated expression, To = {0, 0, lo}, 
So = {0,0, Z }, and Z = — |tanh4^ is a static spin polarization due to the spin splitting 

CO. 

The matrix V can be simplified notably if one takes into account the symmetry properties 
of the Fourier transformation of correlation functions [Eq. (fl31 ]. 72 



(17) 



At the same time, only the symmetrical part T = {{(T^ + r iyM )/2|| is relevant to spin relax- 
ation. This is because the asymmetrical part Ar (T = T + Ar) contributes to Eq. ([161) as 
an additional effective magnetic field that cannot reduce the length of vector AS = S — Sq: 



dAS/dt 



AS dAS 
AS ' dt 



AS 

AS 



[(To - Auo) x AS] - 41 • TAS 



AS 



AS 

AS 



rAS; Aoo x = —AT lJZ , etc. 



Subsequently, the contribution of Ar to spin relaxation is ignored. We also neglect in the 
calculation the second term on the right-hand side of Eq. ( fbil) (proportional to iP) since it 
plays an insignificant role as well.— il&Zi The final result of the symbolic computation for the 
matrix T reads 
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(18) 



Izz + n (lxx + Ixy) -\{lyz ~ Izx) ~ \ly Z 

(ixz ~ lyz) ~ \llz -f ilyz ~ Izx) ~ \ly Z n(^ xx + 7^) + 271%,, J 

1 + e-^)/2, 7 = b>M + J^{u))/2, %„ = i[j^{u) - j Vfi {u)}/2, and 



where n = n(uo) 
7^ = 7^(0)^ 

If the correlation functions possess some supplementary symmetry, Eq. ( TlBl) allows further 
simplification. Assuming 7 M y(u;) = 7^(0;), one can find 



1 ' zz 'yy 
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which readily transforms to the relaxation matrix obtained in Ref. 
In other case with 7 M1/ (u;) = —7^(0;) for \i ^ v, Eq. (fl8|) can be put in a simpler form 
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76 in the limit uj — > 0. 
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(20) 
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Note that Eqs. (fT6|) and (|T8j) are, generally speaking, derived for the case of anisotropic 
medium and arbitrary electron spin splitting u. Equation ffl8l) can be reduced to the common 
Bloch equations if one puts T-f 1 = Tin [y xx {uj) + 7 W (^)], T^ 1 = vr7 zz (0) + T 1 _1 /2 and omits 
the off-diagonal components of T. However this reduction does not always happen, as it 
was shown in Ref. |6J. The spectral properties of Eq. (1171) with respect to the definition of 
n{uj) reveal that the relaxation matrix given in Eq. (Tl8l) is an even function of electron spin 
splitting (i.e., T(uj) = T(— u)), as it also follows from the time inversion symmetry. 77 

The kinetic equations described above do not assume any specific properties of the heat 
bath (Hl) nor the interaction mechanism (Hsl), and form the backbone of our theoretical 
approach. In the following, we show how they can be applied to the problem of QD spin 
decoherence, particularly the elastic processes described by 7^ 2 (0) in Eq. ffTgj) corresponding 
to "anharmonic vibration" and "local field steps" in Table I. 



IV. ANHARMONIC VIBRATION 



There are several manifestations of the elastic relaxation process termed "anharmonic 
vibration". To avoid the difficulties in determining adequate mechanisms of spin-phonon 
interaction and their interference, we consider this problem semi-phenomenologically by ap- 
plying the parameters of the spin-deformation Hamiltonian,— which can be extracted from 
the electron paramagnetic resonance (EPR) experiments on shallow donors in strained crys- 
tals or known dependencies of electron g-tensor shift in a low-symmetry potential. These 
parameters are mainly associated with a direct modulation of spin-orbit interaction as dis- 
cussed above (see Table II). We also take into account the "HFI modulation" mechanism 
when the spin-deformation Hamiltonian is not applicable. A brief preliminary account of 
this process was given earlier.— 



A. Qualitative analysis 

Our analysis is based on the representation of spin-phonon interaction in terms of fluctu- 
ating effective magnetic field Q [in units of energy, see Eq. (jH])] acting on the electron spin 
s. This field is assumed to be composed of additive contributions f2 p from each phonon 
p = {q, k\ with a wave vector q and polarization x, i.e., f2 = S p f2 p . For the moment, let 
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us focus on a single phonon contribution. Then, in the frame of reference rotating with a 
Zeeman frequency, the electron spin performs precession around the small Q p , which oscil- 
lates with the phonon frequency uo p . No alteration in the electron spin phase occurs due to 
such a harmonic perturbation with a possible exception of spin phase shift A<ft acquired at 
the initial period of interaction < t < 2tt/uj p due to a random phonon phase 9 p .— 

A different situation can be realized when the phonon harmonic oscillation is interrupted 
and resumes at a series of instant times tu and hi {i = 1,2,...), respectively, as shown 
in Fig. 3. The reason of such phonon fluctuations can be lattice anharmonicity, phonon 
scattering at the impurities or lattice defects, etc. These irregular phonon perturbations 
affect the electron spin precession resulting in the phase shift A0j at each interval of time 
hi — hi- Subsequently, the net effect of spin phase change <p p (t) due to a phonon mode p 
can be expressed as <f> p (t) = SjA^j, (hi < t). 

Note that for a large number of small changes A0j, their total effect can be described by 
a diffusion equation. Its solution leads to an exponential decay of electron spin phase with a 
relaxation rate T^ 1 = | (A0|) r^ 1 , where r p is the mean time between sequential instants tu 
(or ^2i) -— To estimate the spin phase change A0j caused by a phonon perturbation during 
hi — hi, h is helpful to recognize that the oscillator does not change the spin phase during 
its full period At p = 2tx/uo p as well as for any n integer periods n2Tr/uj p . Hence, A0« can be 
approximated as a spin rotation Q p At p in an effective field Q p independently of the duration 
hi — hi- With the mean value (A0j) on the order of fl p /uj p , one can expect T' 1 ~ T~ l Vt^/ojp 
for the phonon mode p and T 2 _1 ~ SpA^r^ 1 ^/^ when the contributions of all phonons 
(with the population factor N p ) are taken into account. 

B. Theoretical model 

The qualitative consideration provided above shows that electron spin phase relaxation 
can be strongly affected by the phonon phase damping of any origin such as phonon decay. 
Since this mechanism does not involve energy exchange, only the longitudinal (with respect 
to the external magnetic field B) component Q z of the effective fluctuating field is relevant to 
our case. Since the random single spin fluctuation associated with each phonon scattering 
is expected to be very small, we investigate the total result of these small fluctuations 
averaged over the time scale At (r c <C At T 2 ). Thus, the time evolution of the mean spin 
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value s can be described in terms of the quantum kinetic equation discussed in the previous 
section [Eqs. (TIE]) and ffT9]) ]. 64 The analysis of T 2 ~ = nj® z is the focus of the investigation 
as it determines the transversal relaxation stemming from the elastic processes. 

The z component of the fluctuating field operator due to the phonon decay can be ex- 
pressed in a form linear in the creation and annihilation operators at and a p of the phonon 
mode p = {q, x} (q and x are the phonon wave vector and its polarization, — p = { — q, x}); 
i.e., 

n, = ( 21 ) 

V 

Q p = at p - a_ p , (22) 

with a matrix element V? of the spin-phonon interaction, whose specific form will be dis- 
cussed later. Skipping the detailed elaboration of the routine calculation of the phonon 
correlation functions,— 1 ^ the final result for the spin relaxation rate associated with the 
elastic "anharmonic vibration" process is obtained as 

T 27o 1 = 5ZlKf| 2 (2A^ + l)^, (23) 

V P 



where r p is the phonon relaxation time (see Ref. |8l| and the references therein) and the 
phonon population factor N p is given as [exp (hw p /kBT) — 1]~ . This expression is in accor- 
dance with the qualitative analysis discussed in previous section with the exception of the 
second term in the factor 2N P + 1 (i.e., "1") that is attributed to the zero-point oscillations. 
It reflects the quantum effects, which could not be described in terms of classical approach 
but may be important at low temperatures. 

To evaluate Eq. (1231) . one needs the knowledge of the relevant phonon mode p (such as 
the dispersion relation uj p and the relaxation time t p ) as well as its interaction with electron 
spins. Taking into account the conditions frequently adopted in QC, we assume that the 
radius of the electronic state ciq (i.e., the QD size) is much larger than the lattice constant. 
This restricts the summation in Eq. ( 1231) to long wavelength phonons as the matrix element 

falls off quickly for q > 7i/a . Accordingly, r p = r^(q), which is a complex function of 
the temperature and phonon frequencies [see Eq. (6.2) in Ref. [82], can be considered in the 
long wavelength limit. Of a number of potential phonon decay mechanisms, only one term 
originating from the boundary scattering survives at sufficiently low temperatures T < T^ s 
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10 K in the case of Ref. 



821 ). Since this mechanism is insensitive to the temperature, 



81 



it is adequate to assume a constant phonon relaxation time t p ~ r p h- The simple relaxation 
time approximation circumvents the difficult problem of the zone-edge phonons and the 
dependence r p = r^(q), which could be very specific for each particular sample. 



C. Effect of acoustic modes 

To evaluate Eq. (1231) (i.e., |VJ| 2 ) for the acoustic phonons, we consider a phenomenological 
Hamiltonian of the spin-lattice interaction. It can be expressed in terms only of symmetry 
and does not depend on the specific mechanism;— 

(24) 

ijkl ijklm 

Here uu is the strain tensor uja = (dAx k /dxi + dAxi/dxk)/2 averaged over the electron 
ground state \g) (i.e., uu = (g \um\ g)), subscripts ijkl refer to the crystallographic axes, and 
the non-zero components of the tensor AijM can be found from the symmetry of the QD.— ^ 
The second term in Eq. ( 124"]) takes into account the effect of strain gradient wum = duki/dx m , 
which will be neglected from consideration since it can be irrelevant to the elastic spin- 
phonon processes.— 

As a specific problem, we consider a lateral [001] QD with a z-directed magnetic field and 
a localized electron with the axial symmetry with respect to the z axis. This case can be 
applied to high-symmetry systems corresponding, for example, to shallow donors, etc. Now 
Eq. (GU) can be reduced to the form H%% = Os with Q z = [(A 33 - A 31 )u^ + A 31 A\ \lB 
and Q x = Q y = 0. 17 A denotes the dilatation A = TiiUu and the Voigt notation is adopted 
(^4.33 = A zzzz , A 3 i = A zzxx , A m = A xyxy ). Then, as Q z is expressed in terms of the creation 
and annihilation operators [Eq. (I2T]) ]. the matrix element Vf is given as 

( h \ 1/2 

V! = i \jpy^) [(^33 - A 3l ) elq z + 8 H>L A 31 q] $ (q) fi B B . (25) 

Here p is the mass density of the crystal, q = ^Jq x + Qy + q z , Vq is the volume of the sample 
structure, e p the polarization vector of the phonon p, and $ (q) = (g\ e* q r \g). 

The spin-lattice relaxation rate in Eq. (|2"3"|) can be calculated by treating the acoustic 
phonon modes based on the isotropic elastic continuum model with the longitudinal and 
transverse sound velocities vl and vt- Then by taking advantage of the phonon relaxation 
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time approximation (in the long wavelength limit and at sufficiently low temperatures) as 
well as the axial symmetry for the localized electron [i.e., $ (q) = $ (x,z) with x = qao/2, 
z = q z /q], one can obtain 

r°° 7; 3 /'Teff \ 

= t££{B) coth ^-x F» (x) dx, (26) 

Jo x=L,T V * V 1 ' 

(a 33 - A 31 )V|i? 2 

F M {x) = J P„(z)(z 2 + () 2 $ 2 (x,z)dz, 

where we introduce l/a,Q and into £(B) and integrant for convenience. In this equation, 
Pl(z) = (z 2 + C) 2 ; Pt{z) = z 2 {l—z 2 ), is the average phonon relaxation rate as mentioned 
above, Tfff = hv^/kBa is the effective temperature, and £ = A 31 / (A 33 — A 31 ) = —1/3 
if one assumes that the strain induced part of the effective (/-tensor Tjij = YtkjAijkiUki is 
characterized by zero trace, i.e., A 33 + 2A 3 i = 0. 



1. Spin decoherence of shallow donor 



Let us evaluate spin relaxation of a shallow donor electron with an effective Bohr radius 
as (=oo) an d $ (x, z) = (1 + x 2 ) 2 . The last factor determines the relevant phonons, which 
must be long wavelength modes with q < 1/ob- Then, the integral in Eq. fl26l) has a simple 
analytical approximation, 



rp-l 




(27) 



T L V L 

where T"l(t) = (16/15ir)T^U,y this result is for the case of B || [001]. 

As an example, we consider a Phosphorus shallow donor in Si with ag — 1.8 nm. The 
phonon relaxation time can be extracted from the low temperature measurements of Si 
thermal resistivity 8 -^ in terms of the theory developed in Refs. 
10~ 8 s). The spin-phonon coupling constants were estimated in the works of Refs. 



T P h = 2.4 x 
85 and 



86 



However, we believe that direct determination of coupling constants by means of EPR 



measurements of Si:P under an app 
experiment was performed in Ref. 



led stress gives more reliable data. A corresponding 
871 . where the constant A 66 = 0.44 was found. Similarly, 
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our estimation^ obtains A33 = 0.31 and A31 = —0.155 that gives T^q = 1.3 x 10 -4 s _1 at 
the magnetic field of 1 T and low temperatures T <C ^l(t) — 10 K. In another important 
case of a Si shallow donor in Alo.4Gao.6As, the data on EPR under a uniaxial stress^ 8 - provide 
rather strong spin-phonon constants of A33 = 19.6 and A31 = —9.8. This gives the estimation 
T 2 ~ 1 = 6.1 x 10~ 2 s _1 and 6.1 x 10~ 4 s _1 for the magnetic fields of 1 T and 0.1 T, respectively, 
at T = 4 K under the assumption that phonon lifetimes are identical in these crystals. 

2. Spin decoherence in a QD 

Corresponding calculations can be performed for an electron localized in a QD of L xy = 
2ao in the lateral width and L z = eL xy in the thickness. Under the condition e < 0.1, an 
approximate formula takes the form 




where the fitting coefficients are bx = 1, &l = v\jv\, vt = 0.33 — 1.27e 2 , dx = 0.35 — 0.395e 2 , 
VL = 0.97 - 28.5e 2 , and d L = 0.40 - 3.76e 2 . 

To examine the importance of spin phase relaxation caused by the phonon decay, the 
result of Eq. (1281) is compared with the admixture spin-flip process^ in a GaAs QD with 
L z = 3 nm and L xy = 25 nm, assuming r p h = 2.4 x 10 -8 s and A33 = 19.6. For the 
relatively strong magnetic field of 1 T and T = 4 K, the present elastic process and the 
spin-flip give T^q ~ 0.1 s _1 and \T^ 1 = = 10 s -1 , respectively, while for B = 0.1 T 
both predict almost the same rate of ~ 10~ 3 s _1 . At lower magnetic fields, the elastic 
spin relaxation prevails. The Si QD gives a qualitatively similar (but smaller) result; its 
numerical estimation will be discussed in Sec. VII. 

D. Role of optical phonons 

The optical phonons possess high energies compared to the typical energy scale (such as 
the thermal energy ksT) and, thus, are generally ignored in the problem of QC based on QD 
electron spins.— However, this is not the case in the elastic process under consideration 
since the uncontrolled variation of spin phase can happen without the presence of thermal 
phonons as shown in Eq. (1231 at N p — > 0. Zero-point optical vibrations, which survive even 
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at very low temperatures T <C fkj op (u op is the optical phonon frequency), contribute to 
spin decoherence as it happens with acoustic modes in the case T <C T^y [Eq. fl28|) ]. The 
potential significance of the optical modes stems from their relatively large contributions 
at small wave vectors q (|VJ| 2 — > constant as g — > 0) 91 as well as the substantially shorter 
phonon lifetime. In contrast, |V^| — > for acoustic phonons in the long wavelength limit. 

For a detailed analysis, consider a QD made of a polar crystal. The longitudinal-optical 
phonons induce the electron potential 



'27ihw op exp(zqr) ^ ^ 
q 

where \j% — 1/noo — l/ K o! K oo and k are the corresponding dielectric constants. As- 
sume that the QD possesses an asymmetrical potential in the growth direction z due to 
an effective electric field F imposed along this axis (e.g., the gate bias). Since F influ- 
ences the electron ^-factor,— its superposition with ip op results in a spin-phonon interaction 
H°g L = u B B(dg/dF)dip op /dz that is linear in the phonon operator. After substitution of 
Eq. (1291 , this expression reduces to the form of Eq. (ED) with 



Finally, the spin relaxation rate is obtained as: 



T 2 J = r- p % p (B), (31) 

hw op v QD 



'op 

UB) = 4,1^. (32) 



where Vqd is the QD volume (= L L^). 

We perform numerical estimation for a GaAs QD of L xy = 50 nm and L w = 5 nm. Here, 
evaluating the QD potential asymmetry a propri provides the most significant challenge. 
Following the calculations performed for Ga .47ln a53 As/InP quantum wells, 60 a typical value 
for ^rj=r seems to be approximately 0.085 nm/mV at F > 10 mV/nm with B directed along 
the growth axis, while for the perpendicular direction <0.035 nm/mV. As the effect of 
optical phonons quenches in the case of a symmetrical potential (J^ = with F = 0), only 
the upper limit of this contribution can be estimated, i.e., £ op (lT) < 7 x 10~ 6 . Taking into 
account that r^ 1 ~ 10 11 s~\-2i one can find T^ 1 < C op B 2 where C op = 7 x 10 5 s^T" 2 . This 
estimation shows that the transversal relaxation via zero-point optical oscillations can be a 
very effective decoherence mechanism in asymmetrical structures. 
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E. Effects of hyperfine interaction and two-phonon processes 



Let us consider the elastic one-phonon process associated with the modulation of hyperfine 
interaction as well as the two-phonon process. Since these processes result in sufficiently 
slow/long spin relaxation, we shall restrict our consideration to a brief discussion. 

The Hamiltonian of hyperfine interaction with the nuclear spin Ij located at site j with 
the position tj takes the form 

R hf = a hf \M r j)\ 2 !i s > ( 33 ) 

3 

where ip g {r) = \g) is an electron wave function. The spin-phonon interaction due to small 
oscillations at can be expressed in terms of the dilatation operator A = (g\ A \g) if the 
mean internuclear distance (~ rij , ni is the nuclear spin concentration) is shorter than 



the wavelength of relevant phonons; i.e., H^ L = fis, where O = nah/Ay 1(1 + ^uj/Vq^. 
Here, the unit vector n is directed along the effective nuclear field defined by Eq. fl33|) and 
Vqd = (I l^( r )| 4 ^ 3r ) ■ Calculation of the phase relaxation rate for the case of a shallow 
donor results in the expression, which is similar to Eq. ( |27|) . 



-i 

ph 




^ = Tll + ^- (34) 



= I(I + l)nja 2 hf 

Similarly, one can find the approximate rate for an electron localized in a QD through an 
analogy with Eq. (1281) . 

(36) 




= 2.7 — 9.8e 2 , and e < 0.1. Numerical estimations provided for 
a donor in Si and GaAs in terms of Eq. flMl) indicate inefficiency of this mechanism with a 
very long relaxation time (about 10 14 s and 10 8 s, respectively). Hence, this mechanism can 
be neglected in most cases. 

_____ /<2\ 

Along with the Hamiltonian given in Eq. (1241) . the spin-two-phonon interaction Hg£ = 
Yl Dij k i mn fi B BiSjUkiu mn can also contribute to spin decoherence [Dijhimn is the spin-two- 
phonon coupling constants). Now the fluctuating effective field takes the form Q a = 
^2 PP >Wa P 'QpQp'i where W^ ,p ' are the matrix elements of Hg^. Hence the Fourier image 
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of the correlation function 7^ (u) [Eq. (fl3|) ] is expressed in terms of phonon correlation 
functions ((Q Pl Q P2 ) (r) Q P3 Q P4 ) W - Its calculation performed in a harmonic approximation 
leads to a simple approximation for the spin phase relaxation rate 



g factor, deformation potential and energy gap). Numerical evaluation of Eq. ( 1371) at low 
temperatures (T = 4 K) predicts a long relaxation time. In the case of GaAs and B = 1 T, 
one can find T2 5 o ~ 3 x 10 5 s that is too long to be of any experimental or practical interest. 

V. PRECESSION FLUCTUATIONS 

This section analyzes the elastic process listed as "Local field steps" in Table I. The 
specific spin-phonon interactions under consideration are the indirect mechanisms resulting 
from "g-tensor fluctuation" and "hyperfine field fluctuation" (see Table II) via the stochastic 
transitions between different electronic states. 

A. Qualitative consideration 

Consider an electron spin s under the influence of a magnetic field directed along the z 
axis whose strength fluctuates in a step- wise fashion in time (see Fig. 4). In this case, the 
projection of electron spin on the z axis s z is conserved and no longitudinal relaxation occurs. 
Nevertheless, the phase of electron spin will change randomly with the Zeeman frequency 
fluctuation <5f2 resulting in a decoherence rate of T 2 _1 ~ 5Vt 2 r c . Here, r c is the correlation 
time of the fluctuation.— As discussed earlier, the origin of SQ can be the phonon-mediated 
stochastic transitions between different orbital states. Since the spin splitting characteristics 
vary with the orbital states, electrons experience uncontrollable changes in the precession 
frequency (thus, the phase) while undergoing spin-independent phonon scattering between 
them. In the case of a QD, this process may play an important role owing to the shallow 
energy levels and the dependence of g-tensor or hyperfine constant on the electronic orbital 
states.— 




(37) 



The parameter D can be estimated as D = 3(g — 2)C 2 / ' E 2 g (g, C, and E g are the electron 
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B. Theoretical model 



We begin the quantitative analysis by defining the Hamiltonian H over the basis functions 
consisting of a few lowest electronic states \n), which are involved due to the phonon- 
assisted transitions. We also consider that the single-electron problem in the QD gives a 
doubly degenerate energy spectrum E n with the eigenstates \n) in the absence of the Zeeman 
energy and the hyperfine interaction. When a magnetic field is applied, it is conveniently 
assumed that the differences in the spin splitting of the electronic states |n) and \n') are 
small with respect to the energy intervals \E n — E n >\. The specific nature and type of the 
QD is unimportant for the analysis. With these conditions, the total Hamiltonian can be 
expressed as 

H = H s + H e + H p h + H e _ p h- (38) 

The first term H s is the spin (or pseudospin) energy Hamiltonian which can be reduced to 
the form H s = Qs in the most general case with Vt z = (f \H S \ f) — (I \H S \ [} as discussed in 



Eq. 

and 



Its projection on the lowest electronic states |n) reads (see, for comparison, Refs. 
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Q = J2\n)^ n ' n '(n'\, (39) 



where Vt n,n ' are the matrix elements of the effective field (in units of energy) taken between 
the \n) and \n') states. The spin-independent (i.e., orbital) electron energies describe the 
Hamiltonian H e 

H e = Y,E n \n)(n\. (40) 

n 

The Hamiltonians of the lattice and electron-phonon interactions have the usual form 

H ph = ^^ p (a+a p + i); (41) 
p 

H e - ph = J2 B % H ("1 « + a ~ P )- ( 42 ) 

p,n,n' 

Here, B v kk , is the matrix element of the electron-phonon interaction that depends on the 
material parameters and the geometry of the QD. The last three terms of Eq. (j38j) constitute 
the Hamiltonian of the dissipative sub-system Hl (= H e +H p h+H e ^ p h) following the notation 
of Eq. m. 
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Thus, the problem of spin relaxation is reduced to calculating the correlation functions of 
the effective field operator with the Hamiltonian Hl of the dissipative subsystem. Clearly, 
the derivation depends strongly on the specific form of Hl, the energy spectrum, and the 
quantity of electron states considered. Keeping this context in mind, we consider the simpler 
problem of electron fluctuations between only two discrete states |n) = \g) or |e) that are the 
ground and the first excited electronic energy states (with an energy separation of Sq). The 
corresponding Zeeman frequencies are denoted [= QP* in Eq. ([39])] and fii e) (= fi|- c ), 
respectively. Such a simplification allows us to easily perform all the necessary calculations 
in an analytical form. In addition, most of the important physics of the mechanism under 
consideration can be obtained in the framework of this two-level model. 

Hereinafter, it is convenient to introduce Pauli matrices a%, 02, cr 3 on the basis |e), \g), 
where according to the definition, ai, 02, 03 are invariant with respect to the coordinate sys- 
tem rotation in contrast to actual spin matrices s. Then, the Hamiltonian of the dissipative 
subsystem takes the form 

H L = H ph + ^5 a 3 + £ pJ B p 0i(a+ + a_ p ). (43) 

The electron spin Hamiltonian H s can now be split into the steady part H$ and fluctuating 
component H$l with 



and 



H s = ±(n e + W)si (44) 

H SL = ^(n^-n^)sa 3 , (45) 

where 1 represents the unity matrix in the basis |e), \g). According to Eq. (jSj), we must find 
V = H SL - (H SL ) = £ls from Eq. gSJ that determines CI = ±(£1^ - Sl^)(a 3 - (<t 3 )) with 
(C3) = — tanh (5o/2A;bT) . The correlation functions in Eq. ( TlBl now takes the form 

H = i(n<f> - n«)(n(f) - ^)un (46) 

UT) = ((<T3(T)-(a a ))(* s -((T S ))) u . (47) 

The function J^(T) can be obtained by using the double-time Green's function G(t, t') = 
(((73 (t) ; (73 (£'))) with Hl- The final expression takes the following form,— 

UT) = ffi 2 T 2 \ V (48) 
r- 1 = 2^\B P \ 2 {Q P (r)Q- p ) 5l , (49) 
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where the difference between 5' = 5q — \ {Vtz —Vlf ) and 5q will be ignored in the subsequent 
consideration. In the harmonic approximation, one can find 

(Qp(t)Q- p ) s , = (2N P + 1)5 (u p - So) , (50) 

with the phonon population factor N p = ^aj~ a p ) for mode p. The parameter r c has the simple 
physical meaning of the correlation time caused by phonon-assisted transitions between the 
\g) and |e) states. 

C. Temperature dependence of phase relaxation 

Actually, Eqs. ([46l) - (!50l describe the problem under consideration in a very general 
form. As a result, an analysis on the temperature dependence can be provided even before 
the details of the phonon-mediated mechanism that induces the effective field fluctuation 
are determined. Specifically, the correlation time r c given in Eq. ( H9l) can be written as 
t<5 tanh (<5 /2£;bT) with the aid of Eq. (1501) . where t$ = [2ttY> p \B p \ 2 5 (u p — 5q)] is the life- 
time of the excited electron state with respect to the transition to the ground state through 
phonon emission in the limit T —>■ 0. Thus, for the uj = component (e.g., j^ z ), the 
temperature dependence of spin relaxation is reduced to 

J (T) = ^-F (^Yj ! F(x) = (1 -tanh 2 x) tanh x. (51) 

Figure 5 shows the numerical evaluation of Jo{T) as a function of T assuming 5q = 1 meV. 
A maximum is observed at a temperature near ksT = 5q with distinctive slopes on each 
side. The left of the peak corresponds to the reduced hopping from the \g) to |e) state that 
decreases the difference ft^ — (Tl^) (or the amplitude of fluctuations). Hence, in the limit 
T < (5 the effective field fluctuations are frozen and our mechanism becomes ineffective 
as F(5o/2fcsT) — > exp(— <5 //cbT). The slow negative slope on the high temperature side 
[i.e., F (5o/2kBT) — > 5 /2A;_bT] arises due to the well-known effect of dynamical fluctuation 
averaging, which becomes more pronounced with an increase in temperature. 

D. Effect of (/-factor fluctuation 

The general theory discussed above is applied to a specific mechanism of phase relaxation, 
which stems from the hopping between the excited and ground states with different g-factors. 
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In the most general case, the reason for such a difference is the g-factor dependence on the 
energy separation between the discrete electronic levels and the nearest spin-orbital split 
electronic band. For technologically significant III-V compounds, where the interaction with 
the valence band edge determines the deviation of electron ^-factor from the free electron 
Lande factor g rs 2.0023,— one can find the amplitude of the fluctuation Ag = (dg/dE g )5o 
with dg/dEg = (g — g)(A so + 2E g )/E g (A so + E g ), where E g is the band gap and A so the spin- 
orbital splitting of the valence band; we also assume the inequality 5o <C E g . In the case of 
a Si QD, E g is a splitting of the A point in the Brillouin zoneM Then, the final equation for 
the phonon-assisted rate of phase relaxation caused by the fluctuations of Zeeman splitting 
is given by 

T -i_ (g.~g) 2 ( S a (A„ + 2E,) \ 2 ( s a \ 

T ^ - ~W~ {e,(a, + e,)) "° tsF {2^) • (52) 

where huj = g^isB . One can see that our mechanism reveals a quadratic dependence 
of T^~z on the applied magnetic field B in contrast to the ~ B 5 dependence found in the 
previous calculations of longitudinal spin-lattice relaxation through the direct (one-phonon) 



pro cesses . 15 i 16 i 17 i 85 i 86 



An estimation of the excited state lifetime t$ can be performed in terms of a deformation 
potential interaction and a model of lateral carrier confinement. The matrix element of 
the corresponding electron-phonon interaction between the \g) and |e) states is provided 



in Ref. 



17| as B q = iC ^yhq/2pv\\V J osc , where J osc = J osc (q) is a form factor, C is the 
deformation potential, and p, v\\ and Vq are the density, longitudinal sound velocity, and 
volume of crystal, respectively. A straightforward calculation of inverse lifetime results in 



the expression 

T& 32n 2 hr< 2 





,2 



(1 - z 2 )e- a(1 ~ z2) dz, (53) 

ipvj ■ 

where qs = 8 /hv\\, a = 5 /2m e v 2 , and m e is the lateral effective mass. 

To show the efficiency of the mechanism under consideration, we calculate the relaxation 
rates in terms of Eqs. fl52|) and (1531) for a GaAs QD with 5q — 1 meV and the magnetic field 
B = 1 T as a function of temperature (Fig. 6, curve 1). A similar calculation is provided 
for a Si QD with go — 9\\ — 0.0131 and go — g± = 0.0141— (Fig. 6, curve 2). A comparison 
between the result of Fig. 6 and the T\ calculation^ 1 ^ shows that the phonon- induced g- 
factor fluctuation via excited states can control the phase relaxation in a Si QD at T > 2 K 
(i.e., T 2 ,z < Ti) in spite of a small g — g (see Sec. VII for a detailed quantitative analysis); 
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in the case of a GaAs QD, this elastic process prevails over the spin-flip transitions at T > 1 
K and B < IT. 

E. Effect of nuclear field fluctuation 

A similar fluctuation arises in the nuclear hyperfine field when the electron undergoes 
uncontrollable transitions between the orbital |e) and \g) states. As shown in Eq. (15^1) . 
the hyperfine interaction is influenced by the electronic envelope functions evaluated at the 
location of the nuclear spins. Hence, the difference in ^/ e (r) and ^(r) leads to a dispersion 
of the hyperfine field whose mean value is given as 5VL n = dhf\J\l(I + l)xrij /Vqd- Here, 
I and ni are the nuclear spin and its concentration in a QD of volume Vqd as discussed 
before; the dimensionless parameter x = Vqd J (|\I/ e (r)| 2 — \^/ g (r)| 2 ) 2 <i 3 r is equal to 9/16n 
within the approximations of Ref. I17J. If we set Vtf 1 — Vt^ = 5Vt n in Eq. (]4"6l) . we readily 
find the following estimation 



This equation shows the independence of spin relaxation (T 2) hf) on the magnetic field when 
the process is induced by the hyperfine interaction. An estimation of Eq. ( 1541) (with an 
appropriate averaging over 69 Ga, 71 Ga and 75 As)^ for a GaAs QD with a typical size of 
L xy = 50 nm, L z = 5 nm and 5q = 1 meV gives T 2y hf = 8 x 10~ 5 s, 2 x 10~ 7 s, and 
1.7 x 10~ 8 s for T = 1 K, 2 K, and 4 K, respectively. A similar estimation performed for a 
Si QD with the same dimensions and temperatures results in T2.hf = 29 s, 9 x 1CT 2 s and 
6.5 x 10~ 3 s, respectively. Clearly, Eq. (|5"4"|) can provide a dominant contribution for spin 
relaxation at low fields as it does not depend on the applied magnetic field. In the case of 
the Si QD, a detailed quantitative analysis is provided in Sec. VII. 

VI. INTERVALLEY TRANSITIONS 

In this section, we consider one particular but important case that is applicable to the 
semiconductors with multiple equivalent energy minima. The basic principle is the same as 
that discussed in the previous section, i.e., via the phonon-mediated stochastic transitions. 
The difference is that the present case takes into account the intervalley transitions between 




(54) 
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the so-called "valley-split" states. Although normally prohibited, these processes become 
possible in the presence of point defects. Of the two potential mechanisms, the hyperfme 
field fluctuation is suppressed due to the equivalence of the involved valleys. However, the 
" g-tensor fluctuation" mechanism maintains its significance as the degeneracy of these states 
is lifted (i.e.," valley split" with different energies) in asymmetrical structures. 

A. Qualitative consideration 

Following the brief discussion given above, this mechanism can occur in semiconductors 
with multiple equivalent energy minima. As a specific example, a QD grown along the 
[001] (i.e., z) direction with six equivalent minima near the X point is considered (i.e., Si- 
like); the cases with different crystallographic symmetries can be analyzed with a similar 
treatment. Subsequently, the valley-orbital structure of the QD assumes the ground state 
formed by two equivalent valleys [001] and [001] with the wave functions V'ooi and V'ooT- Any 
center- asymmetrical potential V(r) (with respect to the z direction) removes the degeneracy 
over the equivalent valleys with a valley-orbital splitting AE V _ V , thus forming the even 
and odd states tp± = (V'ooi ± V'oolVv^- In the case of a triangular potential V(r), an 
estimation predicts the splitting of AE V _ V = l v ^ v dV (r) / dz with l v _ v m 0.5 A in Si.— Clearly, 
one can see that a moderately asymmetrical confinement potential can result in an energy 
splitting sufficient large to induce the "g-tensor fluctuation" under an appropriate lifetime t$. 
However, the phonon induced transitions between the valley-orbital states with a different 
parity are normally forbidden or sufficiently suppressed.— 

This restriction is lifted when one considers a structure composed of solid solutions. For 
example, the Ge atoms in a Sii^Ge^. QD with a; <C 1 (that offers a promising candidate 
for solid-state QC 96 ) can be treated as the point defects in the Si lattice. The phonon 
scattering with point defects due to the lattice anharmonicity results in the appearance of 
high-frequency harmonics that can transfer electrons between the equivalent valleys [001] 
and [001] or between the valley-orbital states with a different parity. Figure 7 presents a 
diagram of the process responsible for these transitions. As schematically illustrated, the 
third-order anharmonicity splits a phonon resonant with AE„_„ into two high-frequency 
virtual phonons. One of them is "frozen" due to the static deformation surrounding the 
defect, while the other phonon induces the electron to undergo an intervalley scattering. 
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The transfer between the states with an energy separation leads to the fluctuation in the 
electron precession frequency as described in the previous section. 



B. Analysis of intervalley transition 

The estimation of spin decoherence rate due to the transitions between the "valley-split" 
states can be performed in terms of Eq. when the lifetime r$ is calculated taking into 
account the valley-orbital structure of a QD. However, direct application of Eq. (150ft to 
Eq. (|49ft would resulting the vanishingly small matrix elements B q with huj q = AE V _ V for 
the case of intervalley transitions. In order to take into consideration the process depicted 
in Fig. 7, one must supplement the Hamiltonian [Eq. (T4T!) ] with terms accounting for the 
phonon scattering at the defects:— 

H3 = '^jF^qi^i-qQqQqi, (55) 

q,q' 

F q , q >, q i- q = l==i ( 56 ) 

ZpV y/UJqLOqt 

where the matrix element of phonon scattering is 

C q , q > = 20 ^(V,/_ 9 a)(e,a)(e,,a)(qa)(q'a)[(q - q')a]. (57) 

a 

Here gc is a constant, a the unit vectors connecting neighboring lattice atoms, and V g a 



Fourier transformation of the static displacement y(r) caused by point defects. Additiona 
a more general expression is applied for the correlation function (see, for example, Ref. 



ly. 



79 



and 



80l ) in place of Eq. (|50|) that is obtained in a harmonic approximation: 



mr)Q _ p) = l (aVp + PW m 



The phonon damping rate T p (u) must be ca 
Eqs. (|55l)-(l571). Following the works of Refs. 



culated in terms of the interaction specified by 
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and 



971 . it is given as 



T p (u) = 87r \Fp,p',p'-p\ 2 [K^ - u p') ~ ^{°° + Up')). (59) 
<?' 

in terms of the Hamiltonian H ph + H 3 [Eqs. ( HIT) and (|55l) ]. 
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For simplicity, we ignore the lattice anisotropy and consider the displacement y(r) in the 
model of isotropic medium:— >2I 

N D 



y(r) = X>n " "° 3 , (60) 

where is a mismatch in the lattice volume between the host and substitutional impurity 
atoms. As it plays the role of deformation charge, this parameter (i.e., Ap) can be found from 
the analysis of electron paramagnetic resonance linewidth.— The Coloumbic deformation 
field resulting from Eq. (IbOj) gives a well-known expression 

V, = 4.^ -a (61) 

Then, substitution of Eqs. (156]) . (jSTjl and flbT]) into Eq. (159]) and follow-up integration over 
the phonon spectrum in the Debye model (i.e., u p = u q ^ = c^q) results in 



r-R, 



27vr ^ AUh 



^>c, q {oj) = -rr-nDT,^ 7 u u q , (62) 
lb p z c' L 

where njj is the concentration of point defects. The coefficient in this equation takes 
into account the contributions of the phonon modes with different polarizations. In the 
approximation of an isotropic medium, one can find S L = Sll+(cl/ct) 2 Slt, where the sums 
Sxx? = \{X]a( e q^)( e q'^)(^)(^)[('? ~ <?)^] 2 } / account for the polarization mixing; 

\ / (47r) 2 

a = a/a; q — q/g; < ... >(4 7T )i means averaging over the directions of q and q'. A numerical 
estimation for the diamond lattice gives Sll = 0.789 and Slt = 0.239; hence, = 1.29 for 
Si. Finally, the Fourier transformation of the correlation function [Eq. (158]) ] is obtained as 

27 A 2 a 2 in 3 

(Qp(t)Q-p) u = T^( 2N p + l)nnV„-F?—. (63) 

for the phonons connecting the equivalent valleys. 

The next step is to evaluate the matrix element of the intervalley electron-phonon inter- 
action [see Eq. (149]) ]. It assumes the form 



where and S n are the deformation potential constants and the form factor $(q) = 
|e iqr | ip±) selects the wave vectors close to the fc-space separation Ak v of the two minima. 
Assuming a QD with the lateral and transversal dimensions of L xy and L z , one can obtain 

*(<!) « «p[-^(<£ + ql) - ^(q z + AA;,) 2 ]. (65) 
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The appropriate Ak v takes the value of ~ 0.3n/aL (at, the lattice constant) considering the 
Umklapp process. Substitution of Eqs. (162]) and fl64l) into Eq. fl58|) results in the contributions 
of two different types. The first is a direct process and, as mentioned earlier, has only a 
minor significance to r c with a characteristic factor ~ exp[— kf-Ak" 2 ] <C 1. On the other 
hand, the second type accounts for the electron interaction with virtual phonons of energy 
uj q AE V _ V . This contribution is calculated with the help of Eqs. (!63|) -( l65|) : 

r c = r.tanh^^, (66) 



,A 2 n E 2 AE. 



3 



T ' 1 = c *' n ^^? (67) 



When deriving Eq. fl66l) . an estimation g c = 24'-/ pc 2 from Ref. 
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is used that defines the 



constant C iv ~ 210. 



C. Spin decoherence due to intervalley transitions 

With the intervalley transition mediated by the point defects as described above, we can 
now calculate the subsequent spin relaxation due to the associated fluctuation in the local 
magnetic field. Evidently, the hyperfine interaction does not take part in the process since the 
equivalent valley-orbital states possess identical electronic envelope functions. Hence, only 
the g-tensor fluctuations due to the valley split energy is considered. As in Sec. V, applying 
Eq. fl66|) to Eq. fl52|) solves the problem under the condition r c <C T 2 . However, it should be 
noted that the formal calculation of T 2 using these equations violates the required inequality 
condition when AE V _ V is sufficiently small. Subsequently, the Markovian kinetic equations 
[Eqs. ( TT8"j) - (f2"0"|) ] are unsuitable. Actually the spin decoherence time can be determinate by 
r c itself in the regime, where the result for T2 y z [see Eq. ( 1521) ] is much shorter than r c , as 
schematically illustrated in Fig. 8. Here, the spin relaxation can be considered in two stages. 
The first is the "pending" period of approximately r c , during which the electron has not 
undergone the intervalley scattering. Obviously, no change occurs in the electron spin phase. 
Once the transition is made, the electron is subject to the effective magnetic field (due to 
the difference in the g tensor) and starts to acquire the phase change. Within a typical 
time on the order of AT 2 ~ h/ \ — £1^)1 (^ r c ), the electron fully loses the information 
regarding the initial state. The intermediate case T2 ~ r c remains beyond the quantitative 
consideration. For simplicity, we apply an interpolation T dc w T 2 + r c . The calculated 
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decoherence time T& c for a Sii_ a Ge x . QD is plotted as a function of AE V _ V (Fig. 9) and 
temperature (Fig. 10). Two cases of Ge composition (x — 0.1 and x = 0.02) are considered 
assuming the QD size of L xy = 50 nm and L z = 5 nm. Clearly, the proposed mechanism can 
significantly increase the spin relaxation rate in the mixed crystals. Even with the smaller 
"defect" concentration of 2 %, T dc can be reduced to the 10 msec range. Particularly, this 
rate is very strongly dependent on the magnitude of the splitting AE V - V (in reference to 
the thermal energy). Hence, the control of the valley splitting will be essential for long 
coherence. For estimation, we assume H = 8.6 eV, p = 2.33 g/cm 3 , cl = 8.43 x 10 5 cm/s, 



7 = 0.56,-Si and C lv « 210. r 
a Ge atom^ in Si from Ref. 



he deformation charge Ad was found to be 1.2 x 10 25 cm 3 for 



98. 



VII. COMPARISON BETWEEN DIFFERENT MECHANISMS 

To gauge the significance of the elastic processes, the relevant transversal spin relaxation 
rates (T^ 1 ) are obtained numerically in a Si QD (with 5 x 50 x 50 nm 3 ) as a function of 
temperature, magnetic field, and the orbital energy separation. In addition to the three 
elastic processes described in this paper (i.e., anharmonic vibration, local field steps due to 
the g-tensor fluctuation and hyperfine field fluctuation), the contribution by the direct spin 
flip (that is the dominant inelastic process in Si) is taken into account for comparison. 17 Since 
Si is not a mixed crystal and the defect density kept low for the QC application, the process 
via the intervalley transition is not considered. Figure 11 shows the calculated T 2 _1 vs. B 
when the energy separation 5q is fixed at 2 meV. For the range of temperatures considered in 
the study, it is apparent from the result that the direct spin-flip dominates at a sufficiently 
high magnetic field (> 1 T) due to its strong dependence (e.g., B A -B 5 ). For low fields, 
the anharmonic vibration (i.e., phonon decay) process tends to be significant particularly 
when the temperature is also low. As the temperature increases, the contributions due to 
the local field steps become prominent. While the g-tensor fluctuation mechanism plays a 
greater role at a stronger B, the field independent nature of the nuclear field makes it stand 
out in the other extreme. The corresponding map of the dominant process/mechanism in 
the B-T parameter space (5q = 2 mev) is given in Fig. 12, clearly signifying the importance 
of elastic relaxation. 

A similar comparison is provided in Fig. 13 as a function of the orbital energy separation 5q 
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at three different temperatures. The magnetic field is set to 1 T. Note that the contribution 
due to the anharnomic vibration process is not shown as its relaxation rate falls below the 
range of presentation. From the figure, it is clear that the relaxation induced by the local field 
steps (via both the g-tensor and hyperfine field fluctuation mechanisms) are very sensitive 
to Sq; this can be readily understood as the process requires a round-trip transitions by 
overcoming So- Subsequently, its contribution is important only for a small 5q (in reference 
to ksT). Since the spin-flip rate is practically independent of So and T in the considered 
parameter range, the cross-over point between these rates monotonically moves to a higher 
S with increasing T. By 4 K, the inelastic process drops out of the picture. Of the two 
elastic mechanisms, the g tensor case seems to be relatively more efficient as So increases. 
This may be due to the fact that the change in the g tensor is directly proportional to the 
orbital energy separation. Figure 14 summarizes the dominant process/mechanism identified 
in the S -T parameter space at B — 1 T. Considering the conditions typical for QC using an 
electron spin in a Si QD (e.g., S < 1 meV, T < 1 K, B > 1 T), it is expected that T 2 over 
10 sec can be achieved. If a SiGe QD is used in place of Si, then the decoherence process 
can become more active depending on the valley-split energy. However, T 2 on the order of 
10 msec or longer is still attainable. The spin relaxation in a GaAs QD is predicted to be 
much faster. 

Table III summarizes the functional dependence of the dominant relaxation processes 
considered in this study. The zero-th power (e.g., T°, B°, Sq) denotes the independent 
nature to the corresponding parameter. The unique set of T-, B- and (^-dependence will 
facilitate the experimental identification of each process. As the thermal energy provides 
the unwanted error/noise, the parameter space relevant to QC is where T is much smaller 
than the other characteristic energies. 

VIII. CONCLUSION 

A group of elastic spin relaxation processes that do not involve the energy exchange be- 
tween the Zeeman subsystem and the thermal reservoir, is studied in semiconductor QDs. 
Contrary to the common perception, these processes can play a dominant role in the electron 
spin decoherence under the certain conditions (magnetic field, temperature, orbital energy 
separation, etc.) as identified in the investigation. Particularly, the calculation results il- 
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lustrate the potential significance of an elastic decoherence mechanism originating from the 
intervalley transitions in semiconductor quantum dots with multiple equivalent energy min- 
ima (e.g., the X valleys in SiGe). The detailed understanding obtained in this work will 
help optimizing the QD based QC systems for long quantum coherence; for example, the 
Si based QDs can offer a spin relaxation time that is much longer than the GaAs counter- 
parts. However, an additional, comprehensive effort is needed for a complete analysis of spin 
decoherence taking into account all the possible contributions and their manifestations. 
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TABLE I: Main contributions to electron spin dephasing (with T|) in a QD spin based QC 
system. The "random local fields" category denotes the processes that lead to phase dif- 
fusion of the multi-qubit system without causing decoherence in the individual qubits (i.e., 
spin). The decoherence or transversal relaxation (T2) has contributions from both elastic and 
inelastic processes for a single spin. The latter involves spin flip transitions leading to the 
energy (or longitudinal) relaxation of the spin states (i.e., T±). The # sign denotes the spin relax- 
ation processes induced by the spin-phonon interactions. HFI symbolizes the hyperfine interaction. 



Dephasing (T*) 


Decoherence (T2) 


Random Local 
Fields 


Inelastic process (Ti) 


Elastic process 


Direct spin-flip* 
Orbach* 
Two-phonon* 
Nuclear spins 


Anharmonic vibration* 
Local field steps* 
Spectral diffusion 
Geometric phase* 


g-factor dispersion 
Dispersion of HFI 
Impurity fields 
Inter-qubit interaction 
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TABLE II: Various mechanisms of the spin-phonon interaction in a manner consistent to both 
elastic and inelastic processes listed in Table I (marked by SOI and HFI denote the spin-orbit 
and hyperfine interactions, respectively. 





Spin-Phonon Interaction Mechanisms 


Spin-orbit interaction 


Hyperfine interaction 


Spin-spin interaction 


Direct 

Admixture 
Indirect 


SOI modulation 

Kronig— Van Vleck 
g-tensor fluctuation 


HFI modulation 

Spin state mixing 
Hyperfine-field fluctuation 


Waller 
Exchange modulation 

Spin state mixing 
Spin-spin fluctuation 
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TABLE III: Functional dependence of the dominant spin-phonon processes on the relevant 
parameters considered in this study. The zero-th power (e.g., T°, B°, 5®) denotes the independent 
nature to the corresponding parameter. The temperature T, magnetic field B, and Zeeman 
frequency u are expressed in units of energy. T* stands for a characteristic temperature, based on 
which the low and high temperature regimes are defined. Eqs refer to the corresponding equations 
given in the text. 



Process 


Eqs 


Temperature 


Magnetic field 


Energy separation 


T -C T* < T 


T <.uj 








Direct spin-flip 






U) 


T 1 


B 5 






Anharmonic vibration 


|28| 




T eff 


T 1 


B 2 


°o 


g-tensor steps 


|52| 


e So/T 


So 


T -i 


B 2 


e So/T 


So 


HF-field steps 


|54| 


e -5o/T 
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B° 


5 2 e- s °/ T 
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FIG. 1: Schematic diagrams of one-phonon spin relaxation process. The vertical lines correspond 
to the localized electrons with orbital n,n' and spin s =|, [ states; the dashed lines depict a 
phonon with the q or q' mode; and the circle represents a direct electron-phonon interaction. 
When the electronic transition due to the interaction with a phonon does not change the orbital 
states (i.e., re = re'), the diagrams describe the direct process of spin relaxation with the absorption 
(left panel) or the emission (right panel) of the phonon resonant with the Zeeman splitting. If n 
and n' corresponds to the ground and excited states respectively, the net effect of two depicted 
processes in tandem represents the Orbach process of spin relaxation. In this case (Orbach), the 
spin flip occurs only in one of the stages (either absorption or emission) while the other is spin 
independent. 
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FIG. 2: Schematic diagrams describing three types of spin flip s — > s' Raman processes under 
the inelastic two-phonon scattering, (a) The spin-two-phonon interaction in the first order results 
in the direct Raman spin-relaxation process, (b) The process of second order on the one-phonon 
interactions assumes the electronic virtue transitions via intermediate states n's"; the interference 
of the processes, which corresponds to the braced diagrams, sometimes leads to the so-called 
Van Vleck cancellation and the distinctive dependence ~ T 9 for the relaxation rate, (c) The 
anharmonic spin relaxation Raman process implies an inelastic phonon scattering due to the third 
order anharmonicity with the creation of a virtual phonon q" , which in turn induces a spin flip via 
the spin-one-phonon interaction. Processes (a) and (c) are characterized by typical temperature 
dependence ~ T 7 for the Raman spin relaxation rate. As explained, the square represents the 
direct spin-two-phonon interaction, while the triangle stands for the inelastic phonon scattering 
due to the third order anharmonicity; other notification are the same as in Fig. 1. 
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FIG. 3: Schematic illustration of a single phonon contribution to the random spin phase shift. 
^inst is the instantaneous effective field of spin precession associated with a phonon. tu, tu, £13, 
. . . , are the random instants of phonon generation, while £21, ^22, ^23, ■ • • , of phonon annihilation. 
The vertical dotted lines depicts the instants of the last full oscillations by the phonon before its 
disappearance. The residual oscillations (on the right of the dotted line) are not compensated and 
produce the net effect (fi p ), which influences the spin as a short pulses with a random amplitude 
at random times. The phase shift due to the influence of Q p is a Marcovian process at the time 
interval <>t c , where the phonon correlation time r c is defined as a mean value of the differences 
tu + i — t\i (thus, r c = t p , the phonon relaxation time). The phase relaxation time is proportional 
to the mean value of the squared phase shift caused by a single pulse (A^ 2 ) and the rate of phonon 



relaxation r p 1 . 
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FIG. 4: Schematic diagram illustrating the process of precession fluctuations and their influence 
on spin phase relaxation. The model assumes sudden changes of the effective field f2 at random 
instances that lead to the uncontrollable phase disturbances. At the scale of time t » r c (r c is 
the mean interval for consecutive quantum leaps in Q), the Markovian process takes place that 
defines the relaxation rate proportional to T 2 _1 ~ (A0 2 ) /r c , where (A0 2 ) is the mean value of 
the squared phase shift acquired between two successive sudden changes of O. Taking into account 
that A<f> « Qt c , one can find the dependence T 2 _1 ~ (A(/) 2 ) r c typical for the process of precession 
fluctuation. 
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FIG. 5: Decoherence factor F(<5o/2A;bT) [i.e., Eq. (pTD) ] for the process of precession fluctuation as 
a function of temperature assuming 5q = 1 meV. 
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FIG. 6: Spin decoherence rate via the g-tensor fluctuation mechanism [i.e., Eq. (|52p ] as a function 
of temperature for 5o = 1 meV and B = 1 T. Curve 1 is for a GaAs QD, while curve 2 considers a 
Si QD with c/o - 9\\ = 0.0131 and g - g± = 0.0141. 
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FIG. 7: Schematic illustration of the phonon-mediated intervalley relaxation process in a crystal 
with inhomogeneous deformations. The [001] and [001] denote the electron states in the two 
equivalent valleys; the circle represents the electron-phonon matrix element, the triangle depicts 
third-order anharmonicity, and the crossed circle is for the point defect. The real phonon q splits 
into two virtual phonons q/ and q— q/ due to the lattice anharmonicity. The latter is accommodated 
by the local deformation of the point defect that can be of short wavelength with |q/| » |q|. 
Consequently, the other virtual phonon can have a momentum q/ large enough to induce the 
intervalley transition. 
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~AT 2 

FIG. 8: Diagram of phase relaxation in the case of inequality T2 <C r c , where T2 is the transversal 
relaxation time in the Markovian kinetic equations. In time of approximately r c , the effective field 
is subject to a sudden large change. This causes relatively fast spin phase diffusion with the time 
scale of AT2. Hence, the duration of the two-step process, « r c , can be associated with the phase 
relaxation time. 



51 



0.0 0.5 1.0 1.5 2.0 

AE v v (meV) 



FIG. 9: Spin decoherence rate via the intervalley transitions as a function of valley-orbital splitting. 
Curve 1 is for Sio.gGeo.i (namely, the defect density of 10 %), while Sio.9sGeo.02 is considered for 
curve 2. 
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FIG. 10: Spin decoherence rate via the intervalley transitions as a function of temperature. Curve 
1 is for Sio.gGeo.i, while Sio.9sGeo.02 is considered for curve 2. 
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FIG. 11: Spin decoherence rate vs. magnetic field for the dominant spin-phonon relaxation processes 
in a typical Si QD (of 5 x 50 x 50 nm 3 ); line 1 — inelastic direct spin-flip process, line 2 — elastic 
process of acoustic-phonon anharmonic vibration, line 3 — elastic process of local field steps due 
to the hyperfine-field fluctuation, line 4 — elastic process of local field steps due to the (7-tensor 
fluctuation. The temperatures are as shown in the figure; the energy splitting 5o between the 
ground and the excited state is fixed at 2 meV. Since a Si QD is assumed without any point 
defects, the process via the intervalley transition is not considered. 
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FIG. 12: Phase diagram of dominant spin relaxation processes in the Si QD in the B-T parameter 
space. The same conditions as in Fig. 11 are assumed (for example, Sq = 2 meV). The unshaded 
region is where the elastic spin-phonon processes dominate. 
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FIG. 13: Spin decoherence rate vs. 5q for the dominant spin-phonon relaxation processes in a 
typical Si QD (of 5 x 50 x 50 nm 3 ); line 1 — inelastic direct spin-flip process, line 3 — elastic 
process of local field steps due to the hyperfine-field fluctuation, line 4 — elastic process of local 
field steps due to the g-tensor fluctuation. The temperatures are as shown in the figure; the 
magnetic field is fixed at 1 T. Since a Si QD is assumed without any point defects, the process via 
the intervalley transition is not considered. 
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FIG. 14: Phase diagram of dominant spin relaxation processes in the Si QD in the 5q-T parameter 
space. The same conditions as in Fig. 13 are assumed (for example, B = 1 T). The unshaded 
region is where the elastic spin-phonon processes dominate. 
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